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ABSTRACT 


A method for determining the complex bulk modulus B and the 
complex mass density p of a fluid contained in a rigid porous solid is 
investigated. The fluid filled solid is contained within a small 
cylindrical cavity capped on each end by identical transducers. This 
method, suitable for an acoustically "soft" transducer, is based upon 
electroacoustic network theory. B and p are obtained from the input 
electrical impedance of each transducer when they are wired in 
parallel "in-phase" and "out-of-phase", respectively. A description of 


the apparatus and preliminary experimental results are presented. 


ili 


TABLE OF CONTENTS 


I, INTRODUCTION.........10. eee eT 
IT... THEORY oo..cccisses ocaicuicts street er octars see cere estes «2 ee S 
A. INTRODUCTION :. ee ie ees coven se aoa at oes ee 3 
B. ACOUSTIC IMPEDANCE IN SAMPLE CAVITY .................. 3 
Ll. Introduction sce ee een oan ee ee 3 
2. Blocked Acousticgimpeqanece. yest 5 
3. ‘Transfer Impedance v4 
4. Long Wavelength Approximation................cccceeee eee eee 8 
5. Identical Trams@ucenstrn:ss-c2-.0 ccc setts sh cccsen esse 10 
6. Fluid. Filled: Porous Sola arise sees seco esee eee 11 
C. TRANSDUCERS eae eee ee ee eee 11 
D. COUPLER-DRIN IER. Sis Ui Mieco seater ce Saree teaser. ees 14 
III. TRANSDUCER CHARACTERIZATION AND COUPLER........ Nye 
DESIGN 
A. BACKGROUND yoyo ee can as os oycusebue ss acceeree 17 
B. TRANSDUCER SELECTION AND PARAMETER. ............... LY 
CHARACTERIZATION 
1. Matched Transducer Identification ........................0.. 17 
2. Source Strempebe c.g ee roe wise ance 19 
3. Transducer Parameter Determination renee costo eas: 8S Ps, I 
a. Electrical Paraimetensris corse nese nest ens ensees concen 21 
b. Driver Stiffmessiece ous ssi caeia tales 23 
c. Moving’ Mass veers reer eee se aes 28 


iv 


Ge Prans@uction Coecmhiciemt o2..c.c. .ccesce ccc ccvcecccccvcccccecs 30 


me ICCM Ale HImINCSISLATICE 2.455 cureeiierios51sessecseieocscecsesss 33 
fe ANAM es MUN CSS fot oss coe ccen ect eam Nec ce ican wikd'deseeeseeeves 33 
Pima OG EG MRCTIOTIMANCE.. cc. re nerss sere ess <ces os vsae sees 36 
C. COUPLER DESIGN AND CONSTRUCTION................0cee00s 39 
Ie OM PIER OMS UMIICHOM 6202.4... 05cceskscesserssgeteibossdeseverecce ses 39 
DP UNIGe NW AMA BEHICTIL SY SUCIN i snced ccc... vecstssssetccssccsseccasees 41 
IV. PRELIMINARY EXPERIMENTS AND RESULTS .................. 44 
Swe NAIVE Oe ON re ettenia do, ose o.sise's se st's ven se cn ccammennsssschcecvesenesece 44 
See OT I ICE) UIBE, voiss css .osasissornecsveneseoscsecssenscccesdeceeeee 44 
il, JDSereraiepe) 2000.0 0c een enn are nen enn +4 
RC SOMAINCS FR CIMIMEMIU 5. ccacsen sce ccwsedaiewe oie svlsiasledeiisenoeesies 44 
3. Cross-Sectional Area Correction Factor...................... 45 
SR OAG Hs ATlGuH KPVETIMICIIL «. 20.550 0s00s00c0cseseaenestcceseeneserese > 48 
ere OIE EOP REMI NUS ois oo vcs si secsecseconscessesseceesecceonsas 50 
1. Repeat of Water Filled Tube Experiment With Half....50 

of Coupler 
SRE SOM AMC GE SIMU 2222 ccc. sonia cck acc cscescecdecstevesseosedes 50 
ee he Ag OINCCIIOM MN ACLON «0.5 c50caceseceeecs css s00sivsieeeseeceoee 51 
Crs ROAGE AACE MCMIMENCS 2..25...0000 60 -ste0ses<eosseseseeeses 51 
2. Experiments With Complete Coupler......................... 51 
Cem Nim ELS Glee selec evs sents ie ecseasvediessscesseicaeeevecesnns ool 
PWM Y A MEMEO Gleaner ian tec tes sie cawde ieee ai eh iiadarsesssacvesess 54 
V. CONCLUSIONS AND RECOMMENDATIONG ...................0005 57 


ACKNOWLEDGEMENTS 


Several persons and institutions have aided me by providing 
support and assistance. In particular, I appreciate the support of the 
Naval Research Lab Underwater Sound Reference Detachment, 
Orlando. I thank Professors Wilson and Baker for their patience and 
guidance. Finally, I treasured the encouragement of my sons, the 


understanding of my wife, and I thank God all this could happen. 


vii 





I. INTRODUCTION 


Geoacoustic models of fluid-saturated sediment contain several 
parameters which depend on the microscopic pore geometry of the 
sediment [Ref.1, Ref. 2, Ref. 3]. These parameters are very difficult to 
determine from the usual measurements of dilatational sound speed 
and attenuation. In particular, the parameters which determine what 
is essentially the complex effective mass density (p) of the fluid 
contained in the pores are among the least well known [Ref. 4, Ref. 5]. 
For this reason it would be very desirable to directly measure the 
complex effective fluid mass density and thereby extract these 
parameters for various typical microscopic pore geometries. To date, 
direct measurements of p have been conducted with helium filled 
sintered bronze samples [Ref. 6]. 

A method is proposed for measuring the complex effective mass 
density of a fluid contained in a rigid porous solid. The fluid-filled 
solid is contained in a cylindrical cavity capped on both ends by 
identical transducers. The method investigated is an impedance 
method based upon electroacoustic network theory. This method may 
also be used to measure the complex bulk modulus (B), and thereby 
obtain the complex speed of sound, using only a small sample. The 
impedance method is suitable for use with acoustically "soft" 
transducers at the ends of the sample chamber. Moving coil, high- 
fidelity tweeters were used in this investigation. B and ry are obtained 


from the input electrical impedance of the transducers, which is 


related to the acoustic impedance presented to each face, when they 
are wired in parallel "in-phase" and "out-of-phase", respectively. If the 
two transducers are identical in electrical, mechanical, and acoustical 
characteristics, data reduction is remarkably simplified. 

The objectives of this research are to (1) identify and characterize 
two identical, acoustically "soft" transducers; (2) develop a compact 
apparatus to measure both B and p; and (3) validate the technique 
using samples for which the desired values are known or easily 


calculated. 


Il. THEORY 


A. INTRODUCTION 

In this chapter the network equations from which B and p can be 
extracted are developed from the input electrical impedance of the 
transducers. This development is presented in three major sections. 
In the first section the network theory associated with the acoustics of 
the chamber is presented. This development is general in form but is 
simplified by use of the long wavelength approximation. The specific 
case of matched transducers is then analyzed, resulting in a simple 
expression for the acoustic impedance presented to the face of each 
transducer. 

In the second section the network theory for the transducers 
selected for the experiment is presented. The relationships between 
the electrical and mechanical parameters are developed. These 
parameters are then modeled in an equivalent electrical circuit. The 
third and final section consolidates the acoustic network of the 
chamber and the electroacoustic network of the drivers. After this 
consolidation is accomplished, circuit theory is then used to derive B 


and p from the measured input electrical impedance. 


B. ACOUSTIC IMPEDANCE IN SAMPLE CAVITY 


1. Initial Development 
The goal of this section is to demonstrate how the two complex 


parameters B and p can be extracted from the acoustic impedance 


presented to the face of the transducers on each end of a fluid-filled 
porous sample in a cavity. A simple diagram to aid in explaining this is 


shown in Figure 1. 
/-—1-—_- 


porous solid 





Figure 1. Simple Chamber. 


In this figure, there are two transducers, or drivers, facing one 
another with the sample in between. The cross-sectional area of the 
drivers are assumed to be equal to the cross-sectional area of the 
chamber. The two port network representing the acoustic properties 


of this chamber is shown in Figure 2. 





Figure 2. Acoustic Two Port Network. 


In Figure 2, all the parameters listed are, in general, complex. A 


complex parameter is represented with a tilda (~) over the symbol. 


The symbols are defined as: 
p = the acoustic pressure at the face of the transducer. 


U = the volume velocity produced by the transducer. 


Zs = the acoustic impedance seen by a transducer when 
the other transducer is held fixed, or blocked. 


T= the transfer impedance between the two transducers. 
The appearance of Ze and T in each branch of the network is a 
consequence of the reciprocal property of the acoustic field and the 
symmetry of the chamber. Expressions for Za and T for can be 
extracted when the appropriate boundary conditions and wave 
functions are applied, as described below. The subscripts 1 and 2 in 
Figure 2 associate the parameter with its respective driver. 
2. Blocked Acoustic Impedance 

The acoustic pressure and volume velocity can be obtained using 
the concept of a velocity potential (©). If the velocity potential is 
known, the acoustic pressure (p) and particle velocity (a) are derived 


using the following relationships [Ref. 7:pp. 104-106] 


u=VO, (ei) 

~ a® : 

P= -Po> ° where Po is bulk density. (222) 
As an example, for a time harmonic, sinusoidal wave in one dimension 

= ~ t-kx 

( u can now be expressed as u ) of the form elie , these 
expressions reduce to 

~ -kx 

eevee (2.3) 

~ (wt-kx) 
and P=joppPoe (2.4) 


In Equations 2.3 and 2.4, ®9 is the magnitude of the velocity 
potential, u is the velocity of the radiating face, mis the angular 
frequency and k is the wave number. 

Recall VER is the acoustic impedance seen by a transducer when 
the other transducer is blocked. Figure 3 represents pictorially this 
condition (for ease of analysis, x=O is referenced to the face of the 
blocked driver and the analysis performed for the left driver). The 


boundary condition at x=0 requires the velocity there to be zero. 





ele 0 
Co Moving Blocked 


Figure 3. Blocked Acoustic Impedance. 
Applying this condition results in the following forms for the velocity 
potential, acoustic pressure and particle speed (the periodic time 


dependence is suppressed): 


® = Mo cos (kx). (2.5) 
p= -jOPgPo cos (kx). (2.6) 
u=-k@®o sin (kx). (2.7) 


At x= -L these become: 


® = Mo cos (kL). (2.8) 
D = -j® pg Po cos (kL). (2.9) 
u=k ®psin (kL). (2.10) 


For the transducer at x=-L, 
Die Poel... (22.819 
U wus 


where S is the cross sectional area of the transducer face. 


Substituting Equations 2.9 and 2.10 into Equation 2.11 gives 


~ _ -J©PoPocos (kL) 





0“ “Sk @osin (kL) - ae 
Using the relationship kc=q@ yields 
~  -jpoc 
Zq= SPOe cot (cL). (2.13) 


This expression for Ze will be simplified later when the long 
wavelength approximation is applied. 
3. Transfer Impedance 
T is a transfer impedance and is the ratio of the pressure on 
the face of a blocked transducer at one end of the chamber to the 
volume velocity of a moving transducer at the other end (again, x=0 is 
referenced to the blocked transducer and analysis performed on the 


left driver). Figure 4 below represents this situation. 


bear 





0 L 
Blocked [_] Moving 


Figure 4. Transfer Impedance. 


The boundary condition at x=O0 requires the velocity there to be 
zero. Applying this condition results in the forms expressed in 
Equations 2.5 through 2.7. Pressure at x=0 and the particle velocity 
at x=L can then be expressed as 

p(0) = -jopo%o , (2.14) 
and u(L) = -k ®ysin (kL) . (2.15) 
By inspection of Figure 2, note U2 = - S u(L), so that 


~ 


F< = PLO) _ _p(O) __ (2.16) 
Uo -S u(L) 


Substituting Equations 2.14 and 2.15 into Equation 2.16 yields 


ae JO PoPo 


- Sk Oo sin (KL) © (241) 


This can further be simplified to the following expression: 


a “JPoc 


iS S sin (KL) ° (2.18) 


The long wavelength approximation will now be applied to Ze and T. 
4. Long Wavelength Approximation 
This experiment will be conducted at frequencies for which the 
wavelength is much larger than the size of the cavity. The sine and 
cosine terms are expanded in Taylor series form, and since the 
product (kL) is small, the terms of fourth order and higher are 


negligible. The sine and cosine functions are expressed as: 





3 

sin(KL) = - a ; (2.19) 
2 

cos(kL) = 1 - us : (2.20) 


Applying these approximations to Equation 2.13 give 











|. (kL 
~ - -]@ 
Zo= — cot (kL)= J ee —— (2.21) 
Sk L l- (kL) 
6 
Using a Taylor expansion on the denominator i.e., 
ay" Z 
ye (2.21a) 
6 6 
the above equation is reduced to Equation 2.22: 
_ (kL)? 
x _JPoc 3 
Zo= —_———. 
0) S kL (2222) 


Recalling w=kc and B (bulk modulus) = pc?, 


a B jopL 
£0= Sn cee < (27-23) 


Similar relationships are applied to the expression for T found in 


Equation 2.18, resulting in: 


B _jmpL 


Sac GS 


We 


(2.24) 


With these expressions, the acoustic impedance presented to 
the face of each transducer can now be determined. The next step is 
to apply these equations to the case when the two transducers are 


identically driven in-phase and out-of-phase. 


5. Identical Transducers 
When matched transducers are used in the cavity and are driven 

identically either in-phase or out-of-phase, a solution for the acoustic 
impedance (Zac) presented to each transducer face is readily 
obtained. -Consider Figure 2. Using this figure, the acoustic pressure 
at the face for each transducer can be expressed as: 

Pi = ZU, + TU2. (2.25) 

Poe =TU, + ZoUs. (2.26) 
If the transducers are identical and driven in phase, 5 = U 2= U) . 
Equations 2.21 and 2.22 become: 

Py = 20a tp, (2.27) 

po = TU + ZoU = p. (2.28) 
Selecting either equation, and solving for the ratio of acoustic pressure 
to volume velocity gives 

Zae= Zoe. (2.29) 
where Ve is the acoustic impedance. Substituting Equations 2.23 and 


2.24 into the above gives 


Zack aa . (2.30) 
This is the same as for a closed cavity of length L/2. In this case, a 
pressure antinode exists at the midplane of the cavity. 

If the same identical drivers are driven out of phase, the same 


approach is taken and the result is: 


Zao = HORII?) (2.31) 


10 


This is the same as for an open cavity of length L/2 (neglecting 
radiation). In this case, a pressure node exists at the midplane of the 
cavity. 
6. Fluid-Filled Porous Solid 

Equations have been derived for the acoustic impedance of a 
chamber filled only with fluid. _ When sound propagates in a fluid- 
filled, rigid porous body, acoustic energy losses occur due to the 
thermal conductivity and shear viscosity of the fluid. Consequently, 
the bulk modulus and effective mass density become complex. 


Equations 2.30 and 2.31 now become: 


~~ 


See 
Za = OES. (2.33) 


The additional parameter, P, is the porosity. The porosity is defined 
as the ratio of the volume occupied by the fluid to the total sample 
volume. It is required since the cross sectional area for fluid 
interaction is reduced by the frame of the porous material. A formal 


derivation of the porous media effects is presented in the Appendix. 


C. TRANSDUCERS 

Since this research involves determining complex mechanical 
quantities from electrical impedance measurements, it is important to 
first identify the impedance of the transducers used. The transducers 
used are antireciprocal, moving coil type. The theory for this type of 


transducer is well understood [Ref. 8]. 
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Equations have been developed to explain the relationships 
between voltage applied to a transducer and the force on its radiating 


face. The equations are: 


= Zep! + (BI) u, (2.3.4) 


<2 


F = -(BI) I + Zo u. (2.35) 


In the above equations: 


V = the voltage across the input of the transducer. 
I = the input current to the transducer. 
F = the force on the radiating surface. 


u = the velocity of the radiating surface. 
VA = the blocked electrical impedance (V/I; u=0). 
Zine = the open circuit mechanical impedance (F /u; [=0). 


(Bl) = the product of the magnetic field and length of coil 
in the gap of the driver = Transduction coefficient. 


The moving coil transducer can further be modeled by a simple 


equivalent electrical circuit (in the mobility analog) given in Figure 5 





[Ref. 7:p. 352]: 
Ro Ly 
ee 
V 7 
ee 


Figure 5. Driver Equivalent Circuit. 


In the figure above, 


Ro = Blocked electrical resistance, 
Lo = Blocked electrical inductance, 


Ry = Electrical equivalent mechanical resistance, 


2 


Lu = Electrical equivalent mechanical inductance, 
Cy = Electrical equivalent mechanical capacitance. 
The electrical equivalents are obtained from the following 


relationships [Ref. 9]: 





2 
meee) ee (2.36) 
RM Rmech 
Z 
_ (BI) 2.37 
Lm = ( ) 
Cu =——>- (2.38) 
(Bl) 


The parameters Rpecph, K and m are the mechanical resistance, 
stiffness and mass of the driver, respectively. From the above 
relationships, if the stiffness, moving mass, mechanical resistance, 
transduction coefficient and blocked electrical impedance are known, 
the electrical impedance of the driver as a function of frequency can 
be predicted. 

Additionally, for the transducers selected, there is an extra 
stiffness existing due to the air in the cavity behind the dome of the 
voice coil. This stiffness should exhibit the characteristics of an 


acoustic cavity [Ref. 10] given by, 


ke ee (239) 
where 


y = ratio of the constant pressure to constant volume 
specific heat capacity for a gas, 
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S = the area of the driving face, 
Po = the ambient pressure, 
V = the cavity volume. 

This added stiffness modifies the equivalent electrical circuit 
presented in Figure 5. The modified circuit is illustrated in Figure 6. 
The new inductance, Le , accounts for the additional dynamic 
stiffness. Note this model does not reflect any acoustic loading on the 
transducer. The effects of acoustic loading on the electrical 
properties of the transducers will be presented in the next section of 


this chapter. 





Figure 6. Modified Driver Equivalent Circuit. 


D. COUPLER-DRIVER SYSTEM 

When the electrical and mechanical characteristics of the two 
drivers are combined with the acoustical characteristics of the 
chamber, the two port system is given in Figure 7. In the figure 
below, the Voge are the electrical equivalent impedances of the radiation 
impedances presented to the transducers. The radiation impedance 


as used here is a mechanical impedance, i.e., a ratio of force to 
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Figure 7. Combined Electrical Circuit. 


velocity. This differs from the acoustic impedance by a factor of S?, 
where S is the moving area of the driver. 

When the drivers are driven identically in parallel "in-phase" or 
"out-of-phase", Figure 7 reduces to the equivalent circuit in Figure 8, 
where I /2 is the current to flow through each driver and V is their 


common voltage. 


oe 


— <I 


ne 2 
soe ‘ 


Figure 8. Combined Electrical Circuit For Identical Drivers. 


Therefore, by measuring the input electrical impedance of the parallel 
combination (V/I), and knowing Vas and ve Ze can be determined, 


and from it, the value of Tae extracted according to: 


(2.40) 
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Zu Zp 


Defining Z 227. Zo , (2.41) 
ZM + ZR 
_ Bye 

then ZR == (2.42) 

eZ 
The acoustic impedance is extracted from aes in the mobility analog 

2.2 

a Bl) S 
by, Zac= (BI) S_ . (2.43) 
ZR 


Once Dae is determined, the complex bulk modulus is extracted 


using Equation 2.32 when the drivers are in-phase: 


B= jo 





5 BS Zee (2.44) 


or the complex mass density using Equation 2.33 when out-of-phase: 


P= yay Zac: (2.45) 

In summary, by measuring the input electrical impedance of two 
identical drivers used as the end caps of a small cylindrical chamber, 
the complex parameters B and p of the fluid in a porous rigid material 
can be determined, provided the characteristics of the drivers are 


known. 
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II. TRANSDUCER CHARACTERIZATION AND COUPLER DESIGN 


A. BACKGROUND 

The experimental portion of this investigation was conducted in 
two major phases. The first phase involved the selection and 
characterization of two matched transducers. Since it was crucial to 
characterize the transducers as accurately as possible, several 
experiments were conducted, and where possible, alternative 
methods used to verify previous results. 

The second major phase involved the design and construction of a 
test chamber (coupler). The coupler was used to hold the selected 
transducers and the test sample. Once this was completed, 
measurements were made to verify the use of the input electrical 
impedance to extract physical parameters. 

In this chapter the experimental apparatus used in each phase is 
described. Accompanying the description is an explanation of how the 


apparatus functions. Finally, sample data are presented. 


B. TRANSDUCER SELECTION AND PARAMETER 
CHARACTERIZATION 
1. Matched Transducer Selection 
Matched transducers can be obtained by either constructing the 
transducers from identical components or by obtaining many samples 
of a manufactured product and testing several units until two matched 
transducers are found. The latter method was selected, primarily due 


to the experience of Hofler [Ref. 11]. Hofler had worked with 


M7 


acoustically "soft" (moving coil) transducers. The type of transducer 
he used would meet the requirements of this experiment. 

The transducers are nominal one-inch diameter domed 
tweeters. They are manufactured by Philips, Model ADO62T8 
[Ref. 12]. The dome is made of a water-resistant polyvinyl material. 
The cover of the driver is easily removed. Once this is done, the voice 
coil and its support ring can be readily detached from the magnet. A 


cutaway view of the driver (minus cover) is provided in Figure 9. 


o 
Surround ome 


Support Ring 


Cavity 








Yoice Coil 
Figure 9. Cutaway View of Philips Driver. 


Several drivers were obtained and the impedance as a function 
of frequency determined using a Hewlett-Packard Impedance 
Analyzer, Model 4194A (HP-4194A) [Ref. 13], as follows. The drivers 
were mounted on a large baffle and placed in an anechoic chamber. 
The resistive and reactive components of the impedance of each 
driver were measured over the frequency range from 100 Hz to 10 
kHz. The impedance curves for each driver were then compared and 


two drivers with similar responses were selected for further testing. 
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Table 3.1 presents the resistance and reactance values of the 
two transducers selected at their respective mechanical resonance 
frequencies. The similarities in resonance values resulted in their 


selection for further experimentation. 





TABLE 3.1. RESONANCE VALUES FOR SELECTED DRIVERS. 


Frequency (Hz) Resistance(ohms) _|Reactance (ohms) 


954.0 
954.0 70.31 






2. Source Strength 
To further confirm the similarity of the two drivers, the source 


strength of each was found. The source strength measurement is best 
explained by the use of Figure 10. 

The HP-4194A provides the current for the driver. In the Gain- 
Phase mode, it is able to measure the voltage difference between a test 
and a reference channel, or display the voltage of each channel. In 
this case, the reference is the voltage associated with the current 
through a 10 ohm resistor on the return side of the driver. The 
HP-4194A was swept through the desired frequency span. The output 
was then printed on an accompanying printer and stored in the 
impedance analyzer for further analysis. The distance was measured 
from speaker to microphone for source strength calculations. 

The calibration standard was a General Radio (GR) microphone. 
The calibration of this microphone was accomplished by use of a 


pistonphone. The calibration was further checked by using a less 


io 


accurate GR Calibrator. The sensitivity level obtained by these two 


methods was then verified against the calibration curves accompanying 


the microphone. 















HP -4194 






gooooo 





OUTOO GO 
Mic Power Supply 


< To Speaker Test 





(Current) 


10.00 ohm 


Microphone 






Speaker 


a 
}—----— Distance a 


Figure 10. Source Strength Determination. 


Sweeps were conducted for the two selected drivers. Printouts 
of the reference and test channels were obtained for each driver, and 
then compared. When the output from the two different drivers were 


compared, they were virtually identical. Based upon this, these 
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transducers were retained for the remaining portions of the 
experiment. 
3. Transducer Parameter Determination 
Five parameters must be known to adequately predict the 
performance of the moving coil transducers. These parameters are 
the blocked electrical impedance (Zenit driver stiffness (k), moving 
mass (m), the transduction coefficient (Bl) and mechanical resistance 
(Ry). 
a. Electrical Parameters 

The easiest electrical parameters to determine are the 
blocked electrical parameters. Two methods were employed to 
accomplish this. In the first method the voice coil was removed from 
the driver. With the coil removed, a digital multimeter was used to 
determine the direct current resistance of the coil. By use of a 
Hewlett-Packard Low Frequency Impedance Analyzer, Model 4192 
(HP-4192), the inductance at several frequencies was measured and 
then averaged. The second method was similar to the first except the 
voice coil remained in the driver. The moving portion was held fast. 
Resistance and inductance were measured over several frequencies 
and averaged. 

The impedance obtained by these two methods was used to 
determine the best value for comparing modeled to measured 
performance. Table 3.2 presents typical data obtained by these 


measurements. 
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The DC resistance did not vary appreciable from the in- 
magnet case to the out-of-magnet case. However, the blocked 
inductance changed by a factor of ten from the in-magnet to out-of- 
magnet case. The out-of-magnet value was selected for use since it 
provided the best fit to the modeled performance of the lumped 
parameters. The in-magnet inductances were difficult to obtain since 


it was a dynamic measurement and blocking the voice coil was not 


easy. 


TABLE 3.2. BLOCKED ELECTRICAL PARAMETERS. 


| Blocked Resistance (ohms) Blocked Inductance (mH) | 
Coil Removed 
Driver 
aS ae 










6.544 0.11 (at 30 Hz) 


l 
Coil Installed 
ee 






Since it was crucial to know the blocked impedance as 
accurately as possible, the diaphragm of a similar driver with identical 
out-of-magnet electrical parameters was rigidly clamped to the 
magnetic structure using epoxy resin. When the epoxy was cured, the 
blocked inductance was measured as a function of frequency. The 
results are presented in Figure 11. The data reveals the blocked 
inductance is frequency dependent. To properly account for this 
dependence, the blocked impedance over the desired frequency range 


will be obtained from the epoxied driver. The frequency dependent 


Pag 


blocked impedance data can then be stored in the data registers of the 
HP-4194A for later use when the acoustic impedance experiments are 


conducted. 


Blocked Inductance 
160 


140 


120 


100 


wwe ee: 


Blocked Inductance (Henries) 


ee (Hz) 





Figure 11. Driver Blocked Inductance with Coil in the Magnetic Structure 


b. Driver Stiffness 
Several methods can be used to determine the stiffness of 


the transducers. The first method attempted was a static method. 
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The motivation for making a static determination is the ease of 
measurement. If a known mass Am is added to the dome and the 


voice coil deflection Ax is measured, the stiffness can be determined 


by: 
Amg = kAx. (3.1) 


In the experiment performed, the voice coil deflection was 
measured with a linear variable differential transformer (LVDT). The 
LVDT operates with alternating current (ac) and has the property that 
the voltage across the secondary is linearly proportional to the 
position of the plunger. If the change in LVDT secondary voltage is 


measured and the LVDT calibrated, Equation 3.1 can be written as: 


dx 
Amg=kAV 7 , (3.2) 


where AV is the output voltage of the LVDT and dx/dV is the 
reciprocal of the slope of the curve (line) of output voltage versus 
displacement. 

The LVDT used was manufactured by the Collins Corporation, 
Model LMS-119V22-01 [Ref. 14]. It was calibrated using a calibration 
stand and a micrometer. The calibration stand held the transformer 
securely in place while the transformer's lightweight plunger was 
allowed to rest atop the micrometer. As the micrometer position was 
changed, the output of the LVDT secondary was recorded. The data 
obtained were plotted as LVDT secondary output voltage versus 
micrometer position (Figure 12). The slope of this line represents 
the change in voltage for a change in position. Table 3.3 lists the 


slope and its relative uncertainty. 
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fit smoothly around the voice coil dome. 


[Ref. 15]. 


LVDT Calibration Curve 
600 


400 


200 


LVDT Output (mV rms) 


0.0 0.2 0.4 0.6 0.8 1.0 
Position (in) 


Figure 12. LVDT Calibration Curve. 


TABLE 3.3. LVDT CALIBRATION. = 


Slope [mV (rms) / in] Relative Uncertainty (%) 


996.2 





A special ring was machined from lightweight aluminum to 


were added to the special ring. A 1000 Hz stabilized sine wave was 
applied to the primary side of the LVDT, whose plunger rested upon 
the top of the voice coil dome. The change in LVDT secondary voltage 
caused by the change in position due to the weight of the brass ring 


was measured using a Stanford Research Lock-in Analyzer, Model 530 
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Ring shaped brass weights 


Table 3.4 lists the data points obtained using this method. 















ane ee 







Pe. ke 


Two methods can be used to reduce the stiffness data. Both 
were utilized to provide a cross-check. The first involves plotting the 
change in LVDT secondary voltage as a function of mass added to the 


driver. The slope of this line is: 


- -gdv 
Slope = —*—_. (3.3) 
Solving for k: 
ee a GY (3.4) 
Slope dx 


The second method of reduction involves calculating k 


directly from each data point: 


_ _Amg dv 
k= AUass (3.5) 


These values for stiffness are then plotted versus the force applied. 
The result should be a nearly horizontal line whose intercept reflects 
the value of stiffness for small displacements. 

Both schemes were employed to analyze the data. Due toa 


slight but observable nonlinearity in the results from the first method, 
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the second method was selected for final data reduction. The 
transducers are less stiff at low displacements and tend to be more 
stiff at high displacements. Since operation of the driver in the 
chamber experiment will be in the low displacement region, it was 
believed method two would provide a more realistic value for stiffness. 

Table 3.5 lists the stiffness and relative uncertainty for each 
method of reduction. The difference due to the slight nonlinearity is 


seen by the difference in values. 





| TABLE 3.5 STIFFNESS OF SELECTED DRIVERS 
es 
Reduction Method Stiffness (N/m) Relative Uncertainty (%) 
OS a ee ie 
eevee jae le es el 
2380 


Two 










Another way to determine stiffness is to extract it from the 
resonance frequency (@=2z7f) and the (presumed known) moving mass: 
k 


m 


Since and m are known, this equation can be rearranged to solve for 
the stiffness, k: 

k = m(2zf)2 . (3.7) 
This was the method finally selected because it eliminated a sizeable 
discrepancy that existed when the modeled electrical impedance of a 


driver was compared with its measured electrical impedance. As will 


PAE 


be discussed in the next section, the moving mass of an auxiliary 
driver was measured. Using this mass and the vacuum resonance 
frequency of each of the selected drivers, the stiffness for each 
selected driver was solved for and is listed in Table 3.6. The relative 
uncertainty is derived from the uncertainties of the moving mass and 


the resonance frequencies. 
TABLE 3.6. STIFFNESS OF SELECTED DRIVERS __ 


nn ee tt 
Stiffness (N/m) Relative Uncertainty (%) 


5089 





c. Moving Mass 
Several methods can be used to determine the moving mass 
of the transducers. The first is the method of added mass. With the 


driver operating in a vacuum, the resonance frequency is given by: 


| 2 (3.8) 
mM 


where k is the stiffness constant and m the moving mass of the driver. 


If a mass m, is added to the voice coil dome, the resonance frequency 


now becomes: 


k 
SS AV) m+m ° (3.9) 


By rearranging this formula, plotting the square of the period as a 





function of mass added, the slope of the line is 4x? / k and the 
intercept is (4x2 m) / k. Dividing the intercept by the slope will yield 


the moving mass of the driver. 
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The method of using added masses to determine the moving 
mass of the driver was abandoned because of inconsistent results. 
This was due apparently in the inability to obtain satisfactorily stiff 
bonding of the added mass at frequencies of interest. 

A second method is very simple. Assuming all the drivers 
are identical in construction, the mass of any given driver should equal 
that of another, within the tolerance of manufacturing. Using this 
assumption, an estimate of the moving mass can be obtained by 
measuring the mass of a diaphragm removed from an extra driver. 

This was accomplished as follows. The diaphragm of an 
extra driver was cut away from the support ring. The surround was 
then removed from the diaphragm. The mass of each was measured 


separately, the results shown in Table 3.7. 


TABLE 3.7 MOVING MASS OF SELECTED DRIVERS. 
a al 


Pee a 
Mass of Dome and Voice Coil(g) Mass of Surround (¢) 


0.3086 0.0286 





An initial estimate of the moving mass was taken as the mass 
of the voice coil and its diaphragm plus one-half the mass of the 
surround. The value was 0.2943 grams, with a relative uncertainty of 
3 per cent based on half the mass of the surround. This value was 
later refined by curve fitting input electrical impedance data to an 


equivalent circuit model, described in Section II.C. 


29 


d. Transduction Coefficient 

The transduction coefficient (Bl) is obtained directly as 
described below. The experimental setup is the same as the stiffness 
experiment. The two experiments are compatible and can be 
conducted together. 

The driver was placed on a table facing up. A mass was 
added to the dome of the voice coil, resulting in a certain 
displacement of the diaphragm. A direct current was applied to the 
driver and adjusted to return the diaphragm to its original rest 
position, as indicated by the output voltage of a LVDT. By knowing the 
mass added, the applied force was calculated from F=mg. For this 
known force, the current required to exactly offset it was measured. 
Using the relationship F=mg=Bli, the Bl product was found. Several 
masses were added and the BI product calculated for each. The Bl 
products were then averaged and the standard deviation calculated. 


Table 3.8 shows the data for Driver 1 obtained this way. 


TABLE 3.8 TRANSDUCTION COEFFICIENT DATA. 
Mass Added (g) Current (ma DC) Transduction Coef. (T-m) 
P| ee ee | ae 
14.98 









[38.400 [nn 84.46 | nO 
84.84 256.9 3.237 
MM 
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A second method was used to check the values for the Bl 
product. Using the equation F = Bli = -kx, the derivative of this 
equation is taken with respect to current. This results in the 


expression: 
Bl = -k —. (3.10) 


The stiffness (k) is obtained from by the procedure listed in Section 
III.B.3.b. 

The term (dx/dI) is determined by applying a known direct 
current to the speaker and measuring the LVDT output. Table 3.9 
lists sample data for Driver 1. 


_ 
LVDT Output (mV) 
307.86 

302.57 

299.47 

297.19 

294.08 


102.82 288.30 
144.66 285.10 


46.42 291.79 
256.11 276.60 





Linear regression was performed on the above data, with the current 
the independent variable, and the LVDT output voltage the dependent 
variable. The slope obtained by this linear regression is dV/dlI. For 
Driver 1, the change in LVDT output with respect to direct current 
applied was 0.05609 mV/ma. Knowing this and the LVDT calibration, 
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a 
di 


al 
<i 


(3.11) 


a. o. 
=| > 


Applying Equation 3.5 to Equation 3.11, the result is: 


Am g dV. (3.12) 
(AV) dI 


Bl= 

The transduction coefficient was calculated for each mass 

added, using the respective change in LVDT output. The output of the 
LVDT as a function of current applied was determined as discussed 
above. The Bl values were averaged and the standard deviation was 
used to calculate relative error. This resulted in the value of Bl listed 
in Table 3.10. The value is slightly higher than the one listed in Table 
3.8. This was a result of the slight nonlinearity in stiffness. The first 
method used essentially a lower value of stiffness (displacement was 
always near zero). The second method incorporated a higher value of 
stiffness since the stiffness used was the slope of the force versus 
displacement for a wide variety of displacements. This slope was 
influenced by the increasing stiffness with increasing displacement, 
hence a higher transduction coefficient (displacement was over the 


range of driver motion). 


TABLE 3.10 TRANSDUCTION COEFFICIENT OF DRIVER 1 


Transduction Coefficient (T-m) Relative Uncertainty (%) 
3.580 





As a further check of Bl, a measurement of the B field in the 


gap was made with a Hall probe. The number of turns were counted 
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with a microscope and the length of wire estimated from this and the 
diameter of the voice coil. The result was a Bl product of 2.8 T-m. 
This aided in confirming the value of the transduction coefficient. 
e. Mechanical Resistance 

Mechanical resistance is obtained from measured values of 
electrical impedance near mechanical resonance. At resonance in 
vacuum, the magnitude of the impedance is approximately Ro + Ry 
(the reactive component, jmLo , was negligible) from the equivalent 
circuit, Figure 5. Using this fact, the mechanical resistance is 


obtained by rewriting Equation 2.36 in the form: 
(BI) 


R — 
(Zercl- Re 


(3.13) 


Table 3.11 lists the resonance frequency, resonance impedance and 


resulting R,, for each of the drivers. 


TABLE 3.11. VACUUM RESONANCE VALUES. 


Res. Freq. (Hz) | Impedance (ohms) | Mech Resistance (N-s/m) 
662.00 63.11 0.1883 
635.00 50.71 0.2172 





f. Cavity Stiffness 
From the cross section (Figure 9) view of the driver, one can 
observe a large cavity behind the dome of the speaker. This cavity 
provides significant additional stiffness to the system when operated 
in air. Since the drivers will be operated at various static pressures, 
the stiffness as a function of static pressure must be either known or 


eliminated. 
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An attempt to eliminate the stiffness of this cavity was tried 
first. A small hole was drilled through the pole piece of an auxiliary 
driver. The input electrical characteristics were measured for several 
different hole diameters. Instead of eliminating the stiffness, the 
driver now had additional loading caused by the air in the hole. This 
resulted in a coupled oscillator response which would only increase 
the complexity of the problem of added stiffness. 

Since the cavity stiffness could not be eliminated, the next 
step was to determine the stiffness as a function of pressure. 
Assuming the cavity stiffness is due only to the gas in the cavity, it 
should be proportional to static pressure (see Equation 2.39). To 
determine the pressure dependence of the cavity stiffness, the driver 
was placed in a pressure chamber and the pressure was changed in 
controlled increments. Resonance frequency was recorded as a 
function of pressure. The resonance frequency dependence on the 


cavity stiffness is shown in Equation 3.14: 


gael) (3.14) 


ae )(m + m,) 
where: 
f = measured resonance frequency. 
k = static stiffness of the driver. 
k(p) = the pressure dependent cavity stiffness. 
m = moving mass of the driver. 
m,; = radiation mass. 


Table 3.12 lists the data obtained from this experiment using Driver 1. 
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| ~~—sS TABLE 3.12. CAVITY STIFFNESS DATA. | 









| _ Pressure (in Hg) | Res. Frequency (Hz) | 
eee |) 2799. || 





12.3 888.00 
a 912.00 
-0.51 1016.99 


The square of the frequency was plotted as a function of 





pressure. The result was a straight line (the radiation mass, m,;, was 
negligible for the drivers at these frequencies). Once a linear 
relationship had been verified, the stiffness as a function pressure was 
solved for at each frequency. This stiffness was then plotted as a 
function of pressure. The result was another straight line, expressing 
cavity stiffness as a function of pressure. Figure 13 is a plot of these 
data. The equation of the line was solved for by regression, resulting 


in: 
k(p) = (-159.4 +0.07350 (Po) = . (S215) 


The relative error for a value from this equation is 3.9%. The 
intercept value for this equation was not forced to be zero. The 
negative value is a result of the random errors associated with this 
portion of the experiment. 

To check the validity of Equation 3.15, the entire volume 
behind the dome of an auxiliary speaker was measured by filling it with 


a known volume of fluid. Using Equation 2.39, the stiffness was 
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calculated for a pressure equal to one atmosphere. The value agreed 
to within one percent of the value given by Equation 3.15. 
g. Modeled Performance 

Once all the electrical and mechanical parameters were 
identified for each driver, the parameters were converted to their 
equivalent electrical impedances using Equations 2.36 through 2.38. 
Utilizing a circuit analysis program, the impedance of the circuit in 
Figure 5 was modeled on a computer. The resulting modeled 


electrical impedance values were then compared to the actual 


Cavity Stiffness as a Function of Pressure 
8000 


5 8 8 


Cavity Stiffness (mech ohms) 
° 


20000 40000 60000 80000 100000 
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Figure 13. Cavity Stiffness Pressure Dependence. 
electrical impedance measurements made earlier. The first model 
incorporated the electrical equivalent for stiffness obtained using 
added mass (Section III.B.3.b), and the electrical equivalent for mass 


extracted from resonance frequency (Equation 3.8). The results of 
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this model are given in Figure 14 (only the magnitude of the 
impedance is shown). 

To resolve this discrepancy in the fit, the moving mass of an 
auxiliary transducer was measured and stiffness obtained from the 
resonance frequency (as described in Sections III.B.3.b and III.B.3.c). 
Figure 15 demonstrates the performance based upon the revised 


model. The fit has improved significantly. 


Modeled vs Measured using Original Data 


Original Model 


Measured 


!Impedance! (ohm) 





Frequency (Hz) 
Figure 14. Modeled Performance With Original Mass Estimate. 


A best fit was attempted using the built-in circuit analysis 
program of the HP-4194A. This further refined the values used in the 
model which will be used to extract the acoustic impedance from the 


input electrical impedance. Figure 16 demonstrates the excellent fit 
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of measured versus modeled. A best fit was performed for each driver. 
The values were then averaged to define a single transducer whose 
characteristics would be used for the "identical transducers". Table 
3.13 lists the specific parameters and their electrical equivalents to 
be used in the circuit model (Figure 6) to extract the value of acoustic 


impedance. 


Modeled vs Measured Using Removed Mass 
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Figure 15. Modeled Performance With Removed Mass Estimate. 
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~ TABLE 3.13. MODEL TRANSDUCER PARAMETERS 


Parameter ae Baie a 
Blocked Resistance 
Blocked Inductance 
Mechanical Resistance 
Driver_Mass Same Sel 


Driver Stiffness oe 269 
Cavity Stiffness es oe 





Modeled vs Measured for Revised Circuit Elements 
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Figure 16. Best Fit Model. 


C. COUPLER DESIGN AND CONSTRUCTION 
1. Coupler Construction 
Figure 17 is a cutaway view of the assembled cylindrical brass 
coupler (for clarity, only some dimensions are shown). The main 


parts of the coupler are the center chamber with sample sleeve and 
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two hermetically sealed driver chambers, one on each end. The 
drivers are secured in their respective chambers with a backing plate. 

The dimensions of the coupler were dictated in part by the 
dimensions of the driver and by the size of tubing available for the 
sample sleeves. The inner diameter of the sample sleeving used had 
to be slightly greater than the driving face of the transducer. Once a 
sleeving size was selected, the sample chamber size was based upon 
the outer diameter of the sleeving. The length of the chamber was 
chosen to avoid finite wavelength effects for frequencies as high as one 


kilohertz. 
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Figure 17. Cutaway View of Complete Coupler. 
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The sample chamber and the driver chambers are joined 
together and kept aligned via alignment pins. The entire coupler can 
be pressurized. The pressure in the sample chamber is the same as 
the pressure in the driver chambers, minimizing the difference in 
pressure across each driver diaphragm. Communication between all 
chambers is via a vent line. The system is also capable of being placed 
under a vacuum to degas the fluid in the sample and remove bubbles 
acquired while filling the sample cavity. O-rings are used to seal the 
chamber against leakage of gas from the driver side and leakage of 
liquid from the sample side. Samples are enclosed in a copper sleeve 
for easy installation and removal. The sleeve is three inches long and 
has an inner diameter of 1.50 inches and an outer diameter of 1.66 
inches. 

2. Fluid Management System 

Figure 18 represents the fluid system used to support coupler 
operation. A description of fluid system operation is given below. 

To fill the chamber, the pressure source (pressurized nitrogen 
gas) is isolated. The chamber and the reservoir are then exposed to a 
vacuum. The chamber and fluid are evacuated to degas as much as 
possible. The fluid reservoir is then isolated from the vacuum on the 
main vent side. The main vent valve is then cracked open. The fluid 
is driven by atmospheric pressure into the sample chamber. 

The vent line on the chamber is constructed of high pressure 
nylon braided tubing. A fluid level is established at the lower portion 


of this tubing, just above the connection to the chamber. The main 
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vent valve is shut, and the reservoir is again exposed to vacuum. The 
level in the chamber vent line may lower out of sight, so the fill 
procedure is repeated until the level remains visible in the tubing. 
Once a fluid level is established, the chamber is isolated from 
the vacuum and reservoir. The valve to the pressure source is opened 
and pressure is slowly brought up to two atmospheres. Pressure is 
held there for 30 minutes. This ensures that any remaining gases are 


forced back into solution. The liquids to be used in this experiment 
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Figure 18. Fluid System. 
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are Fluorinert liquids. These fluids are manufactured by the 3M 
Company [Ref. 16]. This type of fluid has been used successfully in 


similar experiments. 
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IV. PRELIMINARY EXPERIMENTS AND RESULTS 


A. INTRODUCTION 

Several preliminary experiments were preformed to validate the 
concept that the acoustic load presented to the face of each driver, 
particularly a moving mass, can be extracted from its input electrical 
impedance. Experiments were performed with a water-filled tube 
open in air, with the empty coupler and with the coupler filled with 
water. One important result of these preliminary experiments was the 
need to include a correction factor for the difference in the cross- 


sectional area of the sample chamber and the driver diaphragm. 


B. WATER-FILLED TUBE 
1. Description of Tube 
Once we had confidence in our ability to characterize the 
driver, an experiment was conducted with a water-filled tube placed 
over the driver to see if the mass of water added could be extracted 
from its input electrical impedance. A driver was fitted with a special 
water-tight cover plate with tube, as shown in Figure 19. The tube 
and cover were machined to the exact diameter of the diaphragm, 
including surround. 
2. Resonance Experiment 
A known volume of water was added to the tube, and the shift in 
resonance frequency was recorded. Table 4.1 lists the resonance 


frequencies for the volume of water added. The lowering of the 
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resonance with mass added is expected. These values will be 


quantitatively analyzed further to determine if the resonance 





Figure 19. Initial Test Apparatus. 


frequency shift is that expected, i.e., does the resonance frequency 
shift predict the correct mass added? 





‘TABLE 4.1. RES. FREQUENCY FOR ADDED WATER. 


Volume Added (ml) Res. Frequency (Hz) 
















—E——E————— ae ee 


ee ee SS 


3. Cross-Sectional Area Correction Factor 
The resonance frequency data were reduced by comparing the 
results of the observed resonance frequency with a predicted value. A 
ratio of Equation 3.8 and 3.9 was established, and from the known 


stiffness and mass of the driver, the added mass was determined. 
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When the calculated masses were compared to the masses added, they 
were half of the expected value for all volumes added. 

The data were analyzed in another manner by plotting the 
period squared versus the mass added and is shown in Figure 20. As 


discussed in Section III.B.3.c., the ratio of intercept to slope yields 


Correction Factor Plot 
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Figure 20. Correction Factor Data 


the mass of the driver. When this was accomplished, the calculated 
moving mass was double the known value. 
This apparent discrepancy in the results can be understood as 


follows. It was assumed the driver acts as a piston with uniform 
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displacement across the width of the tube for a specific current and 
frequency. In reality, the edges of the speaker are fixed so this 
assumption is crude at best. 

A better approximation of the driver displacement profile is 


that shown in Figure 21. The displacement is assumed to increase 


L 


—— = 


a 


Figure 21. Assumed Diaphragm Displacement. 


linearly across the surround and then stay constant across the dome. 
Using this profile, the ratio of the average velocity of the fluid above 
the driver was compared to the velocity of the voice coil, which is 
attached to the dome. The kinetic energy of the water added was set 
equal to the kinetic energy imparted by the voice coil. Referred to the 


voice coil, the effective mass of water added is given by: 





Me = 2 ’ (4.2) 
Vve 


where KE is the kinetic energy of the energy of the water, and vy¢ is 
the velocity of the voice coil. The ratio of me to Mwater depends only 
on the diameters of the tube and diaphragm and the diaphragm 
velocity profile. We call this ratio the "area correction factor" since it 
is essentially the square of the ratio of the tube cross-sectional area to 


the cross-sectional area of the dome. The resulting ratio for the 
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velocity profile in Figure 21 is: 


Me 
= 0.71. (4.3) 
water 


m 

In other words, provided the displacement profile is the one 
assumed, the effective mass of the water presented to the transducer 
is only 0.71 times the actual mass of water added. Since the actual 
profile of the transducer displacement is not really known, the 
observed value of 0.5 obtained from the resonance experiments can be 
considered a "calibration" value of Me/Mwater for the tube. A similar 
calibration will be needed for the drivers in the coupler. 

4. Broad Band Experiment 

Once the area correction factor for the tube was obtained, 
experiments were conducted in an attempt to extract the radiation 
impedance (za ) presented to the diver over as broad a frequency 
range as possible (hopefully 100-1000 Hz). Referring to Figure 8, the 
procedure is to obtain the electrical input impedance of the driver and 
the tube with no added fluid (with single driver operation, the current 
is no longer halved). Assuming negligible air loading, this can be 
considered as the series impedance of Ze + Ton . The blocked 
electrical impedance Ge) is subtracted from the input electrical 
impedance and the results are stored in the HP-4194A. The result is 


the electrical equivalent mechanical impedance, a 


Fluid is now added to the tube and the procedure repeated. 
The result is the parallel combination of age and Ze . Since Za is 
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known, Zr can be solved for (Section II.D.) using register 
manipulation on the HP-4194A. Zi is then obtained from: 


Ze (4.4) 
ZR 


The real part of ve is the radiation resistance, while the imaginary part 


is the radiation reactance. 


For the case of fluid in the tube above the voice coil, the 


magnitude of the radiation reactance is equal to wMe where @ is the 
angular frequency and Me is the effective mass added. The HP-4194A 
has the ability to perform the scalar multiplication required to extract 
the added mass. This mass value can be compared to the mass added 
to evaluate the effectiveness of detecting the added mass using the 
acoustically "soft" transducers selected. 

Using the area correction factor mentioned earlier, the values of 
measured mass versus actual mass agreed to within five per cent of 
each other over 80 per cent of the frequency span between 100 to 
1000 Hz. Over the remaining 20 per cent of the span differences 
larger than five per cent can be attributed to the tube resonances that 
occur in the air in the unfilled portion of the tube. These should be 


eliminated when the chamber is tested full of liquid. 
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C. COUPLER EXPERIMENTS 


1. Repeat of Water Filled Tube Experiments With Half of 
Coupler 


Since the area of the sample chamber is greater than the area of 
the tube used earlier, a new area correction factor for the drivers 
installed in the coupler must be obtained. 

a. Resonance Shift 

The initial experiment conducted with the chamber was to 
establish an area correction factor for each driver, as described in 
Sections IV.B.2 and IV.B.3. Only one driver is tested at a time, with 
the driver and the chamber placed in the vertical position. Water is 
added to the open end of the chamber. Several different volumes of 
water are added and the resonance frequency recorded for each. The 
frequencies and added volumes are listed in Table 4.2. With 50 ml of 
water added, the water level was high enough to reach the vent holes 
in the sample sleeving. Double this volume is a good estimate of the 


chamber volume with both drivers attached. 


[TABLE 4.2. CHAMBER CALIBRATION 
ee | 
{ __—10)__ en oo | 
SR ee.) 


146.2 
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b. Area Correction Factor 
The data in Table 4.2 were reduced for the effective mass 
added as described in Section IV.B.3. An area correction factor was 
determined for each driver based upon the reduced data. The average 
of the two area correction factors, to be used in further testing, is 
0.3108. The area correction factors for the two different drivers 
agreed to within four per cent, reconfirming the similarity of the two 
drivers. 
c. Broad Band Experiments 
Once an area correction factor was determined, broad band 
measurements of input electrical impedance were made for added 
water volumes of 10, 20, and 50 ml. Applying the procedure 
mentioned in Section IV.B.4, the data were reduced, resulting in 
extracted radiation mass as a function of frequency. Figure 22 is an 
example of the results obtained for 10 ml of added water. For all 
three added volumes, the extracted values of mass versus added 
frequency were within 10 per cent of the actual value. Over most of 
the frequency band, the data was within five per cent of the correct 
mass. 
2. Experiments With Complete Coupler 
a. Air Filled 
Prior to conducting a test of the chamber filled with water, it 
was first operated with only air in the chamber. Both modes of 
operation were tested: the in-phase "“push-push" mode, and the out- 


of-phase "push-pull" mode. The resulting plots of impedance 
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magnitude versus frequency are shown for the two cases in Figures 23 


and 24, respectively. 
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Figure 22. Mass Extracted from Impedance for 10 ml of Added Water. 


An analysis of the data reveals resonances of the air in the 
chamber. The push-push data showed resonances at frequencies for 
which the chamber length was an integral multiple of the wavelength. 
The push-pull mode exhibited resonances when the chamber length 


was an odd half-integer multiple of the wavelength. Also, these figures 
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demonstrate the drivers do indeed appear electrically as one driver 


with half the impedance, when operated in parallel. Since the results 


Push-Push Mode With Air in Chamber 


lImpedance! (ohms) 





Frequency (Hz) 
Figure 23. Push-Push Operation for Air Only. 
were favorable for the air-only case, the final experiment performed 


was to fill the chamber with degassed water and from the measured 


input electrical impedance of the drivers wired in parallel out-of- 


phase, extract the mass of the water. 
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Push-Pull Mode With Air in Chamber 


lImpedancel! (ohms) 





Frequency (Hz) 


Figure 24 Push-Pull Operation for Air Only. 


b. Water Filled 
With the chamber filled with water, the HP-4194A was used 
to measure the input electrical impedance (as a function of frequency) 


of the drivers wired in parallel out-of-phase. The acoustic radiation 


impedance Za was obtained as described in Section IV.B.4. Fes was 
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further reduced to obtain the radiation mass of the water as described 
in Section II.D. 

The results are shown in Figure 25. The expected value of 
radiation mass was 50 grams. The agreement is poor. At about 200 
Hz, there is evidence of an unexpected resonance in the system. A 


preliminary analysis of the system leads us to believe this is a 


"Push-Pull"' Radiation Mass 


Radiation Mass (grams) 





Frequency (Hz) 


Figure 25. Radiation Mass of Fluid Filled Coupler. 


55 


resonance of the mass of fluid in the vent hole and the spring of the 
air in the tubing above it. It is expected that this unwanted resonance 
can be eliminated by installing a valve in the vent hole which can be 


closed after filling with fluid. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A method has been described by which the complex bulk modulus 
B and complex effective mass density ry for a fluid contained in a rigid 
porous solid can be extracted from measured input electrical 
impedance of a pair of identical transducers wired in parallel in-phase 
and out-of-phase. 

A pair of moving coil transducers was selected and their 
electromechanical performance was characterized. 

A cylindrical brass coupler chamber was designed and built. 

Several preliminary experiments were conducted to verify the 
ability to extract the density from measured input electrical 
impedance. The results for the most part were accurate to within five 
per cent of expected values. The only measurements that exceeded 
this limit were those taken during preliminary parallel driver 
operation for the water filled chamber. Complete characterization of 
the fluid filled chamber still remains to be done. 

Other areas for future research include adding a valve to the vent 
line, repeating the water filled experiment, then interpret the results 
for both push-pull and push-push operation. Once system operation 
has been verified with water, B and ry of simple rigid porous solids can 
be obtained, providing additional validation of the theory prior to the 


testing of the more complex porous solids. 
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APPENDIX 


FLUID FILLED POROUS MEDIA EFFECTS 


This appendix presents a discussion of the effect of placing a fluid 
in a rigid porous solid on its density, bulk modulus and speed of 
sound. Derivations for a simple bulk fluid are followed by those for a 
fluid-filled porous solid. The frame of the porous solid is assumed to 
be rigid. The losses in the fluid saturated porous solid are assumed to 
be associated with the shear viscosity of the fluid. 


The continuity equation for a bulk fluid is [Ref. 17:pp. 106-108]: 


2 + vipa (A. 1) 


For a fluid contained in a rigid porous solid, this formula becomes 


Z <0? rani) a0) (A.2) 





where, 
(p) = the spatial average of the microscopic fluid mass 
density. 


Ge = the spatial average of the microscopic fluid 
momentum density. 


Assuming a one dimensional, time harmonic solution to these 


equations, they become in linearized form: 


jop + poses O, (A.3) 


jo(p) + po=0. (A.4) 
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Now consider Newton's second law for a fluid [Ref. 17:pp. 111-115): 


pats i Bale = -Vp + nVxVxu : (A.5) 
Assuming a one dimensional, time harmonic solution and neglecting 
viscosity for the bulk fluid, this becomes in linearized form, 


~_ 9 
jopou=-<. (A.6) 


The corresponding equation for a fluid-saturated porous solid, 


including viscosity is: 


JOP eg = - ioe Rone (A.7) 


where 
Pere = the effective fluid mass density, 
Cu» = the spatial average of the microscopic fluid velocity, 
(p> = the spatial average of the microscopic pressure, 


Rgow = the flow resistance. 


Perf and Raow are operationally defined by: 


2 
=P ett = kinetic energy density within the fluid, and 


2 
<u» Rpow = rate of decay of fluid kinetic energy density 
due to shear viscosity. 

In general, perp 2p. Also, for steady flow, Ragow = Rpc and 

Rpg = ——!-—5 , (A.7a) 

(pore size} 
where 7 is the fluid shear viscosity. 
It is convenient to combine the velocity terms in Equation A.7 and 


define the complex fluid mass density, p, by: 


99 


jou p p= te, (A. 8) 


Hence 


Row 


16 (A.9) 


p= D ef( aaa 


To develop a wave equation, first the time derivative of the 
continuity equation is subtracted from the spatial derivative of 
Newton's second law. Applying this to the bulk fluid case, 

o~ 
2 
= ae (A.10) 


eo 


For the porous solid, the corresponding equation is: 


(A.11) 


~ o/~ 
2 () 
a © (a9 SP? 
Po ax 
The next step is to momentarily examine the bulk modulus. The 


definition of bulk modulus is [Ref. 18]: 


0 
B= po|[> 


aah (A.12) 








For the fluid in the porous solid the corresponding definition is: 


po P2) 








(A.13) 


Substituting these expressions into Equations A.10 and A.11, the 


following results: 


Q 
noi 


| 


p 


(A. 14) 


Q 
N 


x 


for the bulk fluid, and 
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ee 
2 
-o ‘ey (A.15) 
P Ox 
for the fluid saturated porous solid. 
Equation A.14 is the usual wave equation for the propagation of 
sound in a bulk fluid, with a wave speed given by: 
2b 
c =—. 
p 


Equation A.15 is the wave equation for the propagation of sound in a 


(A.16) 


fluid-saturated, rigid porous solid, with a wave speed given by: 


z 


c = (A.17) 


DU OM 


Hence the speed of sound in a fluid saturated, rigid porous solid is a 
complex quantity which is obtained from the complex bulk modulus 
and complex fluid mass densities in exactly the same way the speed of 
sound in a bulk fluid is obtained from the bulk modulus and fluid mass 


density. 
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